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A UNIFIED THEORY OF PLASTIC BUCKLING OF COLUMNS AND PLATES 

By Elbridqe Z. Stowell 


SUMMARY 

On the basis of modem plasticity considerations, a unified 
theory of plastic buckling applicable to both columns and plates 
has been developed. For uniform compression, the theory, 
shows that long columns which bend without appreciable 
twisting require the tangent modulus and that long flanges 
which twist without appreciable bending require the secant 
modulus. Structures that both bend and twist when they 
buckle require a modulus which is a combination of the secant 
modulus and the tangent modulus. 

INTRODUCTION 

Tlxe calculation of the critical compressive stress of 
columns and of structures made up of plates is an important 
problem in aircraft design. Formulas for the critical com- 
pressive stress have been worked out for a multitude of cases 
of both columns and plates, but these formulas are accurate 
only if the buckling takes place, within the elastic range of 
the material. In present-day designs, most buckling occurs 
above the elastic range. The usual method of handling this 
problem is to retain all the formulas derived for the elastic 
case, but to tiy to discover an effective, or reduced, modiilus 
of elasticity which will give the correct result when inserted 
into these formulas. 

Column buckling was the first structural problem to be 
studied in the plastic range. In the latter part of the nine- 
teenth century, Engesser proposed use of the tangent 
modulus as the reduced modulus for columns. At almost 
the same time, in the belief that the column would be 
strengthened by unloading on the convex side^ Consid&ce 
suggested that the effective modulus should lie between the 
tangent modulus and Young’s modulus. This concept was 
subsequently refined by Engesser and by Von Klrmdn 
(reference 1) and led to what is generally knovTi as the 
“double modulus.” 

Experiments have shown, however, that the Von Kdrm^n 
double modulus gives values that are too high for the column 
strength (reference 2) and that the correct modulus is 
probably the tangent modulus. Shanley (reference 3) has 
stated the situation compactly as foUow's: “If the tangent 
modulus is used directly in the Euler formula, the resulting 
critical load is somewhat lower than that given by the re- 
duced modulus theory. This simpler formula, originally 
proposed by Engesser, is now widely used by engineers, 
since it gives values that agree very weU with test data.” 
Further careful tests by Shanley (reference 4) and also by 
Langley structiues research laboratory- have shown that tjde 
unloading on one side of the column, postulated by 


Von K4rm4n, does not occur at buckling and that the correct 
modulus for columns is actually the tangent modulus. This 
conclusion also has theoretical justification (references 3 
and 4). 

In the case of local or plate buckling, the reduced modulus 
is appreciabl 3 ' higher than the tangent modulus. Tests 
of the local buckling stress of aircraft-section columns have 
been made b\' Gerard (reference 5), who has suggested the 
use of the secant modulus for tins type of buckling. Exten- 
sive tests in the Langlej’ structures research laboratory on 
similar aircraft sections made and reported over a period 
of several years and summarized in reference 6 have also 
shown that the reduced modulus for plates is in the vieinitj- 
of the secant modulus. In particular, tests of long aluminum- 
alloy cruciform-s^tion columns, designed to buckle bj’ 
twisting without appreciable bending, have been made 
in a manner similar to that described for the aircraft-section 
cohunns in reference 6. The results have Shown that the 
reduced modulus for pure twisting is verj" close to the 
secant modulus. 

The present paper constitutes, a theoretical investigation 
of the buckling of plates be^’ond the elastic range, which 
includes columns as a limiting case. Such an investigation 
requires a knowledge of the relations between the stress 
and strain components beyond the elastic range. These 
relations have not as j^et been conclusively determined. 
A recent paper bj" Handehnan and Prager (reference 7) 
based on one possible set of stress-strain relations led to. 
results for the buckling of hinged flanges in sharp disagree- 
ment with test results obtained at the Langley structures 
research laboratory. Another set of stress-strain relations 
is generally accepted by the Russian investigators and has 
been applied bj- Ilyushin (reference 8) to the stress conditions 
in thin plates. 'These results form the foundation of the 
present paper, which assiunes that in plates as weU as in 
columns unloading during the early stages of buckling does 
not occur. On this basis, a unified theory of plastic buckling 
applicable to both column and local buckling has been 
developed. The results are presented in the following section. 

RESULTS AND CONCLUSIONS 

Hyushin (reference 8) has treated the stability of plates 
stressed above the elastic limit with consideration of the 
three possible zones that might result from buckling: 
(1) a purely elastic zone, (2) a zone in which part of the 
material is in the elastic and part is in the plastic state — 
the “elasto-plastic” zone, and (3) a purelj" plastic zone in 
which all of the plate is stressed bejmnd the elastic limit. 
AU three zones may exist simultaneously if the plate is not 
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entirely in the plastic state before buckling or if the buckling 
is allowed to proceed beyond the initial stages. 

If, however, the plate is uniformly loaded before buckling 
so that all parts of it are initially at the same point in the 
plastic range and ifj in addition, buckling and increase, in 
load are assumed to progress simultaneously, then the plate 
may be expected to remain in the purely plastic state in the 
early stages of buckling. This second assumption is in 
agreement with the corresponding condition that apparently 
holds for columns (reference 4). 

Upon the assumption that the plate remains in the purely 
plastic state during buckling, Ilyushin’s general relations for 
this state have been used to derive the differential equation 


of equilibrium of the plate under combined loads. Since 
critical stresses are generally simpler to compute from energy 
expressions than from a differential equation, the correspond- 
ing energy expressions were also found. Those derivations 
are given in appendix A, together with applications to com- 
pressive buckling of various types of plates. A comparison 
with Ilyushin’s treatment of the plastic-buckling problem is 
given in appendix B. 

The remits of most interest in the present analysis are 
given in the following table as values of a quantity ij, the 
number by which the critical stress computed for the elastic 
case must be multiplied to give the critical stress for the 
plastic case. 


Structure 

t ■ ■ 

n 

Curve 
(See fig. 1) 

Long flange, one un- 
loaded edge simply 
supported 

E _ 

A 

Long flange, one un- 
loaded edge clamped 

%s(0.330 + 0.670yi+|^) 

B 

LohYplates, both un- 
loaded edges simply 
supported 

/l , 1 /l,3E,an\ 

T^'\2'^2V4‘^4 

C 

Long plate, both un- 
loaded edges 
clamped 

%= (0.362+ 0.648^1+11^) 

D 

Short plate loaded asa 
column 

1 Etan 

4 E ■'■4 E 

E 

Square plate loaded as 
a column 

0.114 ^-hO.886 ^-2 

F 

Long column 

. E~ ' - 

— — ■■■ — ^ ^ . . . ■_ _ - . 

G 


These values of rj are plotted as curves A to G in figure 1 for 
extruded 24S-T aluminum alloy for which the compressive 
yield stress was 46 ksi. Similar curves for a could readily 
be prepared for any other material having a known stress- 
strain relationship. 

The values of 17 given in the table were obtained by 
dividing the critical stress of the structure in the plastic 
region by the critical stress that would be obtained on the 
assumption, of perfect.. elasticity. Since Poisson’s ratio has 
been taken as one-half in both computations, errors. from 
this cause will ordinarily be present in both critical stresses. 
Most of these errors will be eliminated, however, in the 
process of division to obtain tj; and, consequently, the values 
of ri given are believed to be nearly correct. 

When plate-buckling stresses in the plastic range are to be 
computed, the experimental value of Poisson’s ratio that 
applies as closely as possible to the stressed material, to- 
gether with the appropriate value of t) from this paper, 
should be used in the plate-budding formula. 


E 

The highest value of 17 which is can be realized only if 

there is negligible longitudinal bending (as with a long hinged 
flange which buckles by twisthig). The low^est value of 17 
E 

which is occurs when the longitudinal bending pre- 
dominates over other types of distortion (as with a long 
column under Euler buckling). The theory implies that a 
change in the stress^train curve caused by prestressing of the 
material would alter the value of 17 in the first case but not 
in the second^ if the buckling stress is higher than the highest 
stress reached during the operation of prestressing. If, on 
the other hand, the buckling stress is lower than the highest 
stress reached during the operation of prestressmg, then 17=1 
for each case. 

liANGLEY -MEMORIAL AERONAUTICAL LABORATORY, 

National Advisory Committee for Aeronautics, 
Langley Fieldv Va., "JuTy 29 , 1947 . 
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Compressive stress, Jtsi 

Figure 1.— Compated coives shovinc varfatloa o( :i with stress (or various structures of 24S-T slumlnum sUo7 lu compresslau. (Curves A to O are drawn for a material with a yield 
stress of 46 ksL Circles are test data from cmelform sections and the scatter band enclosing curve A shows the limits of variation of the cruciform-section properties from 46 ksl. 
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APPENDIX A 

THEORETICAL DERIVATIONS 


Definitions.' — ^The intensities of stress and strain are 
defined in reference S', respectively, as . 

a-j='Vo-,‘'+o-„*— o-jO-ir-j-Sr* (1) 

VsV ^**+‘»*"^***‘''*'^ ■ ^ 

where 

ffr stress in the r-direction .. . 

ci strain in the x-direction 

ffj, stress in the j/-direction . . 

ty strain in the ^/-direction 

T shear stress 

T shear strain 


in which ei and «2 are middle-surface strain variations and 
et is the middle-surface shear-strain variation, Xi and xt 
are the changes in curvature and x* is the change in twist, 
and 2 is the distance out from the middle surface of the plate. 

The corresponding variations 5Sx, SSy, and St in Sx, .Sy, 
and T must be computed. From equations (3), 

Sx=E„cex 

therefore .. 

5iS'i=£’sec5«x+eiS 


According to the fundamental hj’pothesis of the theory of 
plasticity, the intensity of stress <r< is a uniquely defined, 
single-valued function of the intensity of strain for any 
given material if at increases in magnitude (loading con- 
dition). If at decreases (unloading condition), the relation 
between <r< and a* becomes linear as in a purely elastic cas& 
In the equations of definition (1) and (2),- the .material is 

taken to be incompressible and Poisson’s ratio=|- The 

stress-strain relations compatible ■with the equations of 
definition (1) and (2) are: Z. . 



(3) 


3t 



Now the variation of the work of the internal forces is 

a-iSet—axSex-\-aySey-\-TSy 

so that 

atSfx+<ryS€y+TSy 
» 

(Ti 

a^xfi+<fyfi+2rft—ziaxXi + ayXa+iTXs) 

Vi 


Substitution of this value of Set in equation (6) gives 
SSx=Ei,^5ex— 

[Vi€l-4-Vj«2-i-2T«3— 2(<rxXl + VyXj-|-2TXj)] 

Let the coordinate of the sm-face for which S«i=0 (the neutral 
sui-face) be z=Zq. The expression for 2 o is obtained by set- 
ting 3e<=0_in equation (6); 

ax,jp^,+ 2re, 
v,Xi + v,Xj+2rXa 


By inti’oduction of this coordinate into the expression for SSx 

and^'asE,.,. 


and by recognition of — as 


det 


These relations imply isotropy of the material. 


Variations of strain and stress. — "When buckling occurs, 
let €r, €y, and y vary slightly from their values before buckling. 
The variations Sex, Sey, and Sy ^vill arise partly from the 
variations of middle-surface strains and. partly from strains 
due to bending; thus. 


ei— sxi" 

Sey= e»— 2 X 2 '■ 
57=2«3— 22 x 3 . 


( 4 ) 


SSx=Egec(fl — SXl)-|- 

(Aieo— E'tan) (ffxXl + V,X2+ 2rXi) 2o) (8) 

a fit 

In a similar way it may be shown that 
SSy=E,te(ei — 3 x 2 )-!- 

; (^»eo— ■®t«i)(v*Xl + ffrX2-t-2TX»)(2— * 0 ) (9) 

a fit 

and 

2 

Saeo (es — 2Xa) + 

(fi'.ec — -Et»n) (v*Xl + V»X2+2rXs) (3~ 2o) . (10) 
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Variations of forces and moments. — ^For the variations of the impressed forces Tx, Ty, and Txy and the moments 
My, and Mxy 

k 




&M, 


SM, 


'x=j\saxdz 

2* 

k 

J 2 

5cr/fz 

“T 

■=r 
=/ 


dr dz 


■z dz 


I SffyZ dz 


bMx 




5tz dz 


on 


whore h is the thickness of the plate. 
From equations (3), (8), and (9), 

SA/^I J (sSx+^ SSy)z dz 


_ ^ e*+- £ 

■Elee + ^ ea^J* ^sdz — Fiev- ^Xi + ^ X2^J* ^^dzi - £'t,n) (o-*Xl + <7,X2+ 2 tXj) J ^ (z* — ZZo)(/z 


‘3 12 






( 12 ) 


where 

Similarly, 


r\/ 

9 


SMy=-^J' ^ ^SjS'^+^ bS^ z dz 

— ®' { [‘-I Ci)’ (> -fc)] >“+5 [' -I (i -fc)] X--I ,T (i-fe) 

k 

1'2 

5A/i^= I ^brzdz 

— T { 0-fe)] X.) (i-l^) } 


(13) 


(14) 


In these expressions, the integrations of bSx, bSy, and 5 t in the plastic region have been taken over the entire thickness of the 
plate, with the implication that no part of the plate is being unloaded. 
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Equation of equilibrium. — If v) (x, y) is the bending deflection of the plate at buckling, and if no external moments are 
applied to the plate^ then the equation of equilibrium of an element of the plate may be ivritten 




bxdy 


b^(SMg) I / b^w , b^w , 


5% \ 
bxbyj 


in which the impressed forces <r*A, tryA, and rh are considered as given (virand are positive for compression). 
w, the changes in curvatures are 

and _ . . . _ . _ . 


The change in twist is 

_ 

^~bxby 


(15) 

In terms of 
- (Ifia) 

(10b) 

(16c) 


When the values of 5Mx, dM^, and SAfi, in equations (12), (13), arid (14), respectively, are differentiated as required by 
equation (15) and substituted in that equation, the general differential equation of equilibrium for a plate in the plastic'statc 
is obtained as follows: 


— [-fl 3v,<r,+ 2r*/ 


b*w 

M 

by*^' 

D'\ 


/ b-w. 


b'^w 


bxbyj 


b’‘w \ 


In the elastic range, equation (17) reduces to the usual form 

vv— ^( 


b% , b^w \ 


where 






(17) 


Energy integrals. — ^Equation (17) is the Euler equation that results from a minimization of the integral 


r (Y—irf rr b^,^ r/b% Y , b’w b*w"l „ 

J J \ 2 t ' \SzV * ba:“ dxby"^ *L\^^/ * 


hai I 
~2 


b^W ^ f^'^W 
bxby * \byv J 


(18) 


which represents the difference between the strain energy in the plate and the work done on the plate by {the external forces. 
The coefi&cients in this integral are : 


In the 'plastic region ... In the elastic region 


«=‘-l(7:J0- 

'E,J 


c,=\ 

a=z^(i-^' 

) 


c,=o 

^ , 3 VjV„+2t^ 

(‘ t:) 

' 

Ci=i 

C4=3^(l-fi5' 

) - 


( 74=0 

‘"•“■-J (?,)’(*- 

■^taa\ 

■®«ec/ 


( 7«=1 
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If there is a restraint of magnitude e along one longitudinal edge of the plate, the strain energy in this restraint itself is taken 
to be 


1 ^ f Vdt£\ “ ‘ 

2 6 J Lv2>2//»->'ol 


dx 


(19) 


if 3/o is the edge coordinate. (See reference 9 for form of expression.) In expression (19), the stiffness D' is assumed to be 
the same as that in equation (12) . If restraints are present along two edges, there will be two terms similar to expression (19) . 
These terms may be added to int^ral (18) as additional strain energy. 

Critical stress in plastic region. — ^If the integral (18), supplemented if necessary by additional terms of the form of ex- 
pression (19), is set equal to zero and the resulting equation solved for ct, the critical-stress intensity in the plastic region 
(o-i)piis 

_D' JJl ^bx^'dxby *bxbyc)y^ 1 ^ L\^/y-yiJ 








( 20 ) 


in which the values of the O’s in the plastic range are used. This expression for the critical-stress intensity may be minimized 
as with the corresponding elastic case. 

If the values of the CPs in the elastic region are used in formula (20), the critical-stress intensity in the elastic region 
(ffj), I is as follows: 

rmnzr: 

^bxj 




( 21 ) 


Expression for ij. — quantity i; is defined as 


(^<)p» 


(22/ 


This quantity is a direct measure of the effectiveness of plasticity in reducing the critical stress of a structure, and its com- 
putation in terms of the constants of the stress-strain curve represents the solution of the problem of plastic buckling. 

Application to plates compressed in the x- direction. --The theory will now be apphed to flat, rectangular plates uni- 
formly compressed iu the x-direction. Values of ij will be computed for the following cases: 

I. Long plates with one free edge (flanges), the other edge being either binged or’ clamped 

II. Long plates with both edges either hinged or clamped 

III. Plates with two free edges (columns) 

"When o-„=T=0, and the plasticity coefidcients reduce to 


^ 1 I ^ 

<^ 1 - 4+4 X.0 


(?!,=<7«=0 

Cj=Cj=l 

The differential equation of equilibrium, equation (17), then becomes 

^ d*w , „ dho I d*w h(Ti ihj} 

Z?' bx* 

and the corresponding energy expression (20) for the critical stress in the plastic range becomes 


(23) 


(vi) 


Pl' 




^ \bx*/ ~^\jbxbyj 

* I 5*W r I 



[dx 

* J 

r/dw\ 

JV£)a:> 

2 

) dx dy 


(24) 
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Case I : Elastically restrained flange _ 

If 2/=0 is the elastically restrained edge of the flange and y=b is the free edge, a deflection surface known to be good in 
the elastic range and presumably satisfactory also beyond this range is (reference 9) 


where 


[(!)+- (!)+- (!)+»■ (!)],*- T 


Ui= -4.963 
U 2 — 9.852 ---- 

flj= -9.778 

and € is the magnitude ot the elastic restraint. Substitution of this expression for w in equation (24) gives 




(Vx)pl — ^2 


06-2^:- 


M 


+ 




where 


1 I Cje , ^ 

3‘^2a,‘^4o3' 


W 


Ci=0.23694 . 
Ca=0.79546 
Cb ^0.89395-.. . 

In order to find the minimum value of ((Xx)pt, 


which gives 


The minimum value of (<rr)pj is therefore 


C4= 0.04286 
<•,=0.56712 
<•8^0.17564 


C7=0.I9736 
C8=— 2.3168 
<^=4.0982 


(?)■ 




(5+¥+'-S^)^‘ 


— _2 
7T 




1 , Cse I Cae^ 

3 ”^ 203 '^ 403* 


irD' 

hVi 


For the elastic case, the same expression is obtained from equation (21)withCi=l andZl' replaced b.yZ?. From equation (22), 
therefore, 




+ T)(i+T+^') 


(25) 


(a) If the edge y=0 is liinged, e=0 and^ from equation (25), 


E 


(26) 


This value, as a function of sti’ess, is plotted as curve A for 
24S-T aluminum alloy in figure 1. The individual points 
represent the NACA tests of the buckling of ..crucifonii- 
section columns for which the. condition e=0 is fulfilled. 

(b) If tire edge y=Q is clamped, e=oo and, from equa- 
tion (25), 




i 0®”^ C7^+2 v^i 

i(c8-icr) 


Vc*c, 
Vc^ 


- ^=%(o.330 + 0.670^|+|^) (27) 

This value of tj is plotted as curve B for 24S-T aluminum 
alloy in figure 1. 
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Case II: Plate elastically restrained along two unloaded edges 

Ify=±^ are the immovable unloaded edges which are 

elastically restrained against rotation by restraints of magni- 
tude e, a satisfactory deflection surface is known to be 
(reference 10) 

t] X 

Substitution of this expression for w in equation (24) gives 
(<rr) 5 ,=[c, 0 4-/i(e) Q) +/ 2 (e)]^ 


where 


and 




/.(e) = 


0.0237eH0.297«-|-^ 
0.0046 l£^-f0.0947(=-l-| 

0.0114e2-f-0.1894e+l 


0.00461e*-|-0.0947«-t-^ 
In order to find the minimum value of 

=0 


^ (a~ i) 


which gives 
The minimum value of (or*)j,j is therefore 


(xHW 




(ffi) j 1 = [2 VCi -\/i (e) +/ ! (e) ] 

For the elastic case, the same expression is obtained from 
equation (21) with <7i=l and D' replaced by D. From 
equation (22), therefore, 


_ -^«ee 2 y Cl -//l (e) 4~/. (e) 


(28) 


(a) If the edgesy=±2 are hinged, e=0,/i(«)=l,/i(€)=2, 
and, from equation (28), 


Case III: Plastic buckling of columns 
For the discussion of the plastic buckling of columns, it is 
convenient to revert to the differential equation (23). The 
plate, when loaded as a column, has two free edges described 
by the conditions 

/'bho , 1 bhD \ 

W 2 2)a;V _ . ^ 

ir— ±2 


\Py^ ' 2 c)3^c)y/ 


s-±-. 


A solution of equation (23) which identically satisfies the first 
condition is 


w=(^2 cos I cosh ^4-p cosh ^ cos cos 


I 


where 




/3= "--y/ j 7 + "\/*+ (t) 

2=^^+KtJ 

1 I 3 ExKa 

^i-4+4 

In order that this solution also satisfy the second condition 
at the free edges, it is required that 


tanh 


i'. 2-i fan _ 




which is the buckling criterion for the plate when loaded as 
a column. Let 






(32) 


.Eiec 1+ V^_.Eaec /I , 1 /l , 3 £'un'\ 
E 2 £ V 2'^ 2^4"^ 4 £'„<■/ 


(29) 


This value of 17 is plotted as curve C for 24S-T aluminum 
alloy in figure 1. 

(b) If the edges y=±.^ are clamped, e=«>, /i(e)=5.1.5, 

/2(e)=2.46, and, from equation (28), 

E,„ 2.46-|-4.52Va 
E' 6.98 

=%(0.353 + 0.647.^/fi|^) (30) 

This value of 17 is plotted as curve D for 24S-T aluminum 
alloy in figure 1. 


where 5* is a quantity to be determined for three individual 
cases. By use of equation (32) 
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and the buckling criterion given in equation (31) becomes 


c?q; 


tank ^ 

a 

2 


or 


O 

tan ^ 

■|3V -^^0 


tank ^ 




(i-Vi-f*)] 


o 

tan I 

~T 

Q 


= 0 


(33) 


From equation (32), ^=(j^ and thus the critical 

stress in the plastic range is . . 

(<r*)ji- 3^Zy 


where 


4' 

h 


^ Vl2' 

The coiTesponding critical stress in the elastic range is 


(<r.).j^ §77 



The reduction factor ij is obtained from formula (22) as 

.=a-£i 


I ^00J_3 


(34) 


(a) In order to investigate the case of short columns, let 
f approach zero. Then, by definition of /3, 


/3-^0 


and 


tan^ 


1 

2 




In addition, 




The buckling criterion given in equation (33) therefore 
reduces to 


tank ^ 


[,+(i-.)(.+v.-?)]:7i=o 


The expression in the brackets approaches ^ as f->0. In 
order to satisfy the buckling criterion, therefore, 


tanh ^ 


a 

2 




which can_be realized only if a is large; that is, if b/l is large. 
For short columns, therefore. 




and, from equation (34), 


1 I 3 

4 ~E~^4t - E~ 


(35) 


This value of )j is plotted as curve E for 24S-T aluminum 
alloy in figure 1. 

(b) For a square plate, ^1, a=xVl + Vl-^ 

)3=ix-\/ 1 — Vl — f®, and the buckling criterion given in 
equation f 33) becomes 

tanli ^^1 + Vl~ f 






tanh (|Vi—7r^) 
which is satisfied by 0.1 5375. From equation (34), 




V=0.114 ^+0.886^ 


(36) 


This value of 77 is plotted as curve F for 24S-T aluminum 
•alloy in figure 1. 

(c) For long columns, « and P become so small that 


tanh ^ tan ^ 

• Of — - — m 1 


a 

2 


2 


and the buckling criterion of equation (33) reduces to 

(i~ f”) (1 - 0 


which is satisfied by From equation (34), 


Etio 


(37) 


This value of 17 agrees with the experimental results of 
references 2 and 4 and is plotted as curve G for 24S-T 
aluminum alloy in figure 1. 



APPENDIX B 

COMPARISON WITH ILYUSHIN’S STABILITY CALCULATIONS 


The basic difference between Iljnshin’s solution of the 
plastic-buckling problem and that given in this paper is that 
Ilyushin considers the plate to unload on one face as it 
buckles. The unloading process results in the creation of 
an elastic-plastic zone in the plate, and different equations 
from those that apply when the plate remains plastic during 
the buckling process are required for this zone. 

The differential equation for the buckling of a rectangtilar 
plate when buckling is accompanied by unloading is given 
by Hyushin as equation (3.43) of reference 8. For simple 
compression in the a;-direction this equation is of the same 
form as equation (23) of the present paper, but with the 
following different constants; D is used instead of D' and 

A=l-Xf*(3-2f) 

E 

is used instead of In the formula for k, from equa- 

-^aec 

tion (3.1) of reference 8, 

.. i-VT=x 

X 

and, from equation (1.22) of reference 8, 

X=l-^- 


When the values of f and X are inserted into the expression 
for k 



E, 


Computation shows that k is always larger than so 

xl/sec 


E, 


that the use of t in place of will result in appreciably 
higher values of than those given in the present paper. 


Since Ilyushin uses the elastic value D, there is no possibility 
of the solution yielding a secant modtilus. Curves A to G 
in figure 1, if computed from Ilyushm’s equation (3.43), 
would start with a horizontal line at unity for curve A 
(Young’s modulus) and end with curve G expressing the 
Edrmdn double modulus which is appreciably higher than 
the tangent modulus of this paper. If D' were substituted 
for D in Ilyushm’s equation (3.43), curve A would then 
represent the secant modulus as it does in the present paper, 
but curve G would still remain the Kirmdn double modulus. 
Therefore, when the unloading of the plate during the 
buckling process is considered, results are obtained which 
are not confirmed by experiment. 
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